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Abstract 

This letter is concerned with the analysis of the six-vertex model 
with domain-wall boundaries in terms of partial differential equations 
(PDEs). The model’s partition function is shown to obey a system 
of PDEs resembling the celebrated Knizhnik-Zamolodchikov equation. 
The analysis of our PDEs naturally produces a family of novel deter¬ 
minant representations for the model’s partition function. 
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1 Introduction 


The formulation of the Quantum Inverse Scattering Method (QISM) [lj2 represented a 
large step towards the understanding of algebraic structures underlying exactly solvable 
models of Statistical Mechanics [3j. In particular, it benefited from many insights gained 
through the study of the six-vertex model originally proposed to explain the ice residual 
entropy [4]. The relevance of two-dimensional vertex models to other fields has gradually 
increased over the years and nowadays we can find ramifications from condensed matter 
physics to gauge theories. For instance, the relation between two-dimensional classical 
vertex models and one-dimensional quantum spin chains [51 is the most notorious one but 
recent applications also include the study of supersymmetric gauge theories |6|. In the 
latter case exactly solvable vertex models are among the main players in the so-called 
Bethe/gauge correspondence describing the supersymmetric vacua of two-dimensional 
JV = 2 gauge theories in terms of Bethe vectors |6|. In addition to that, they also play a 
prominent role in the program initiated in 17] enabling integrability based techniques for 
the computation of anomalous dimensions of certain gauge invariant operators in = 4 
super Yang-Mills theory. 

Roughly speaking, the formulation of a vertex model involves three main ingredients: 
graphs on a lattice, statistical weights for graph configurations and boundary conditions. 
Although those ingredients are arbitrary a priori, integrability in the sense of Baxter 
|3,[5] imposes restrictions on all of them. The statistical weights are then required to 
satisfy the Yang-Baxter equation while boundary conditions are usually constrained by 
compatibility conditions [8j,|9]. The latter often consist of algebraic equations but some 
types of integrable boundary conditions are not characterized in such a way. This is 
precisely the case for domain-wall boundary conditions introduced in [l0| as a building 
block of scalar products of Bethe vectors. 

The six-vertex model with domain-wall boundaries has also found several important 
applications ranging from enumerative combinatorics to the study of gauge theories. For 
instance, in 11 Kuperberg demonstrated that the partition function of the six-vertex 


model with domain-wall boundaries exhibits a combinatorial line counting the number of 
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Alternating Sign Matrices 12 . In the gauge theory front, the chiral partition function of 


a two-dimensional Yang-Mills theory with gauge group SU (. N ) was shown to correspond 
to the partition function of the six-vertex model with domain-wall boundaries in the 
ferroelectric regime 13 . The literature devoted to this particular variant of the six- 
vertex model is quite extensive and we refer the reader to 10,14-18 and references 


therein for a detailed account. In particular, its partition function can be expressed as 
determinants 14,19 or alternatively as a multiple contour integral 20-22 . As far as the 


integral representation is concerned, it is worth remarking that it was obtained in 21,22 


as a solution of certain functional equations governing the model’s partition function. 
Those functional relations have their roots in the Yang-Baxter algebra and, interestingly, 
they also give rise to a hierarchy of higher-order PDEs satisfied by the six-vertex model 
partition function 


23,24 


The analysis of this variant of the six-vertex model in terms of PDEs is the main 
purpose of this paper. However, the PDEs we shall consider here are not the ones de¬ 
rived in 24 ; although they originate from the very same functional relation obtained 
in ~~ ~ 


22 


One important feature of these new PDEs is that they consist of linear first-order 

Moreover, 


equations resembling the celebrated Knizhnik-Zamolodchikov equation 25 


we shall demonstrate here how one can simply read off the solution from the struc¬ 
ture of the equations. Through this analysis we obtain a whole family of determinant 
representations as a result of Cramer’s method for solving linear systems. 


Outline. This paper is organized as follows. In Section [2] we describe certain functional 
equations satisfied by the partition function of the six-vertex model with domain-wall 
boundaries. Those equations have their roots in the Yang-Baxter algebra and in Section[2] 
we also discuss some of their properties which will be required throughout this paper. In 
Section [3] we describe a procedure producing first-order PDEs from the aforementioned 
functional equations. Through this approach we obtain PDEs with structure similar to 
that of Knizhnik-Zamolodchikov equations. The resolution of our PDEs is then discussed 
in Section [4] while concluding remarks are left for Section [5j 


2 Functional relations 


Several methods have been formulated along the years aiming to evaluate partition func¬ 
tions of lattice models and the Kramers-Wannier transfer matrix technique 26,27 plays 
a distinguished role. Within the Kramers-Wannier approach, the computation of the 
model’s partition function is converted into an eigenvalue problem for the associated 
transfer matrix. Although the latter might be regarded as an auxiliary tool for this 
computation, it gained a more fundamental status with the advent of Baxter’s concept 
of commuting transfer matrices jsj. For instance, in models possessing commuting trans¬ 
fer matrices there usually exist functional equations characterizing the transfer matrices 
eigenvalues. The largest eigenvalue is then responsible for the model’s free-energy in the 
thermodynamical limit. However, it is worth remarking that a more direct description 
of partition functions in terms of functional equations has also been put forward in 


28 


The six-vertex model with domain-wall boundaries is a singular example where the 
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computation of the model’s partition function can not be readily recast as an eigenvalue 
problem along the lines of Kramers-Wannier technique. Nevertheless, this particular 
model still admits a description in terms of a spectral problem as shown in 29 . As 


far as this variant of the six-vertex model is concerned, the possibility of describing its 
partition function in terms of functional equations was firstly discussed in 30 . The func¬ 
tional equation presented in 


30 


is, however, restricted to values of the model anisotropy 
parameter satisfying a root-of-unity condition. For arbitrary values of the anisotropy 
parameter an alternative functional equation was subsequently derived in 31,22 . The 


latter originates from the Yang-Baxter algebra and it will be the starting point of our 
present analysis. 

Functional equation. Let Z be the partition function of the six-vertex model with 
domain-wall boundaries on a L x L lattice 


10 . It is a symmetric function Z: 


—>• 


depending on L spectral variables A* £ C. Strictly speaking Z also depends on a set of L 
variables /q £ C, usually refereed to as inhomogeneity parameters, but here they will be 
fixed. Now write X l,j := {A*, \ i < k < j} and define the sets X 1 ^ := AU J \{Afc}. Using 
the algebraic-functional method formulated in 22,32 , one can show that the partition 


function Z satisfies the following functional equation 


T M * Z X L ) = o. 


( 2 . 1 ) 


4 = 0 


with coefficients 

( L 


M,, := 


II 6(A 0 - i-ij) - H a(A 0 - i-ij) H T7v~—v4 

7=i 7=i 7=i ax ~ 


c(Aj — Ao) 
KK ~ Ao) 


H a(Aj fij) H 


3 = 1 


Q* (A j A i ) 
K A, - A*) 


for i = 0 

otherwise . 


( 2 . 2 ) 


The functions a, b and c in (2.2) correspond to the statistical weights of the symmetric 
six-vertex model and here we shall restrict our discussion to the rational model. In that 
case we have a(A) := A + r], 6(A) := A and c(A) := r/. We can readily see that a = 6 + c 
which corresponds to a quotient of the six-vertex model algebraic curve a 2 + 6 2 — c 2 = Aab 
with A = 2. The parameter r] is usually refereed to as semi-classical parameter and here 
it will also be fixed. In this way c(A) can be regarded as a constant. 

As a matter of fact Eq. (2.1) does not consist of a single functional equation for 
the partition function Z. It comprises a total of L + 1 equations uncovered through 
permutations Aq £>■ Aj. This feature has its roots in the following properties: 


i. Eq. (2.1) runs over L + 1 variables whilst Z depends only on L variables; 


ii. Analytic solutions are symmetric as demonstrated in 21,22 ; 
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in. Although (2.1) is invariant under permutations A* -H- A j for 1 < i,j < L, the per¬ 
mutation A 0 A i produces a new equation with the same structure but modified 
coefficients. 


The above remarks pave the way for extending (2.1) to the following set of functional 
relations, namely 


£ M<"> Z(X°- L ) = 0 


n — 1, 2,..., L 


(2.3) 


i=0 


with coefficients 


r ( n ) ._ 


1 Ao-H-An 

for i = 0 


M 0|a o ^ 

for i — n 

(2.4) 

M-i 1 Ao-H-A n 

otherwise 



In what follows we shall describe a mechanism rendering a system of first-order partial 


differential equations similar to the classical Knizhnik-Zamolodchikov equation from (2.1) 


and (2.3). 


3 Partial differential equations 


Several analogies between (2.1) and the theory of Knizhnik-Zamolodchikov equations 
have already been pointed out in 24 . However, the different nature of the aforemen¬ 
tioned equations seems to prevent a direct identification at first sight. Although the 
existence of a precise matching is not clear at the moment here we intend to shed some 
light onto possible connections by unveiling a system of first-order PDEs underlying Eqs. 


(2.1) and (2.3). 


In order to proceed, we shall firstly restrict our attention to Eq. (2.1). In Section [2] 
we have already pointed out the distinguished role played by the variable A 0 and we can 


readily see from (2.2) that the coefficients M 0 and Mk exhibit a simple pole when Aq = A 


Despite the existence of such poles the limit Aq \k in (2.1) is still well defined and by 


taking this limit we are left with a first-order PDE satisfied by the partition function Z. 
For later convenience we also introduce the following extra definition. 


„±i „±i 


r-ilT 


be the ring of formal 


Definition 1 (Evaluation map). Let A n := 

Laurent series in n variables over C and let & n denote the symmetric group on n letters. 
In addition to that define A n := 1 , • • • os the subset A n C A n formed by 

symmetric Laurent series. Next set i,j G {1, 2,..., n} and the evaluation map £: A n -A 
A n _i is defined by the plethystic substitution 


( £ijf) (xi,x 2 , ...,x n )\= f(x i,x 2 , ...,x r 


(3.1) 


for f e A n such that the RHS of (3.1) is well defined. 
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Theorem 1. The partition function Z satisfies the system of PDEs 


diZ — f 'y ^ Ojj + hi j Z , 


3 = 1 


(3.2) 


where <9j := and i G {1, 2,..., L}. The terms f2jj and hi are in their turn defined by 

n _ c(Aj — Aj) a(Aj — Aj) -p-|- a(Aj — p£ b(Xk — A*) a(Afc — Aj) c 

13 A7A V\T7\ rrllTTA—ATT 11 ^ 


a(Aj Aj) &(A* Aj) ^ a(Aj /^fc) o(A k A*) &(A^ Aj) 

k^=i,j 


j rr ^(Aj — /ifc) -rr &(Afc — Aj) ^ c(Aj — pfi) c(Afc — Aj) c(Afc — Aj) 

11 TTxTTTTT 11 TTxTATTT + 2^ TTTTTTTTT + 2^ “ ~~ ■ 


^ ^ o(Aj hk) ®(Afc Aj) | a(Aj /^fc) | b(^Xk Aj) a{Xk Aj) 

k^i k^i 


(3.3) 


Proof. Set Aq = Aj + a in (2.1) and expand it in power series in a. By doing so, and taking 


into account that Z is a symmetric function, we are left with a Taylor series containing 


only positive powers of a. Then (3.2) is obtained by setting a = 0. Alternatively, one 
can also take directly the limit Aq —> Aj of (2.1) using L’Hopital’s rule. □ 


The structure of (3.2) resembles that of a classical Knizhnik-Zamolodchikov equation 


and we refer the reader to 33 


34 for more details on such equations and their gener¬ 


alizations. Moreover, Eq. (2.1) is not the only functional equation originated from the 


algebraic-functional framework as discussed in Section [2} It turns out that we can asso¬ 
ciate a system of first-order PDEs to each one of the equations (2.3). Those equations 


exhibit the same structure of (3.2) and they read as follows. 


Theorem 2. For i,n € {1,2, ..., L} the partition function Z satisfies the following 
family of partial differential equations, 


d,z= ( J 2 T”> + h X) Z. 


3 = 1 


(3.4) 


The terms and h^ in (3.f) are explicitly given by 


n 


(n) «(Aj — A j) a(\ n — Aj) t r a{\j — Pk) -r-r b(Xk — Aj) t r a(Afc — A j) 

l] b(Xj — Aj) a(A n — Aj) ^ a(Aj — pfi) ^ a(Xk — A,) b(Xk — Aj) 

k^i k^j 

L L 


£ 


13 


«£? := c-n«- i (A,-mn 6(Al ' A,) 


k =1 


k =1 

k^i,n 




3 ± n 
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h(X n Aj) h(X n hk) T 


k=l 


a 2 (Xi- X n ) 
b(Xi X n ) 


J"Jfl(A n /^fc) J J 


k=i 


o(Afc X n ) 
6(Afc — A n ) 

k^i,n 


, (n) _ c(Aj — ^i) c(Aj — Aj) c(Aj — Hj) a(Xi — X n ) 

1 ^ 6(Aj - Ai) a(Aj - Ai) j-( a(Xi - fij) ~~ b(Xi - A n ) 

iT- 

Proof. Same proof as for Theorem [lj 


- 1 


£. 


(3.5) 


□ 


Remark 1. Eq. (3.4) fori = n coincides with Eq. (3.2). Thus we have a total of L 
systems of PDEs and not L + 1 as one might have expected. We shall restrict ourselves 


only to the system (3-4) through our analysis in order to avoid over counting. 


4 Solution 


The resolution of the functional equation (2.1) has been already described in 22 
the method devised in 


In fact, 


22 seems to apply to several functional equations with structure 


resembling (2.1). In this way we have obtained multiple contour integral representations 


for a variety of partition functions and scalar products associated to six-vertex models 
through this approach |21 22,35-37|. Despite the sound results, the applicability of this 
method still seems to rely on special properties of the desired solution which might not 


exist for certain models. For instance, the method devised in 21,22 depends strongly 


on the existence of the so called special zeroes which, in case they exist, might still be 


difficult to locate. This is precisely the point where the formulation (3.4) seems to be 


more suitable than (2.1). In what follows we shall demonstrate how the solution of (3.4) 


naturally emerges from its structure without requiring strong properties of the partition 
function Z. 

Theorem [2] gives us a total of L systems of first-order PDEs and we can readily notice 


the presence of L — 1 terms of the form Z in the RHS of (|3.4t) . Our equations can 


TJ 


be fortunately solved for those terms and in order to see that we write 


VWZ =: < 4 % 


(Z) 


(4.1) 


in such a way that can be directly read off from (3.5) 


Next consider the subset of (3.4) formed by n G {1, 2,..., i — 1, i + 1,..., L}. We can 
conveniently rewrite the associated L — 1 system of equations as 


5><”’4(Z) = (cd. - ;4”») z 


3 = 1 


(4.2) 


in order to make the dependence with the L — 1 terms £ij(Z ) more apparent. In this 
way we can readily apply Cramer’s rule and by doing so we find 


£u(Z) = 


det (H, 




det(Wi) 


diZ 


det(R, 




det(Wi) 


(4.3) 
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with matrices W % , and H t j defined as 


and 


(i) 


,(!) 


m := 


+1 

^ii—1 

^ii+1 

■■■ +L 

, ,h-i) 
co iX 

, ,6-i) 

U ii -1 

6-i) 

w ii+l 

, ,6-i) 

• • • W iL 

. ,6+i) 
io a 

, ,6+1) 
U ii -1 

, ,6+i) 

^ii+l 

, ,6+b 

• • • W iL 

K t ; (L) 

\ w il 

,+) 
U ii -1 

,+) 
W ii+1 

,+) ! 

• • • w iL / 


w <?> 


Hi, - 


i 1 


CO. 


co, 


6 - 1 ) 

il 

6 + 1 ) 

il 


CO, 


( 1 ) 

*3-1 


CJ. 


( 1 ) 

*3+1 


ca 


(<-i) 

27 — 1 
6 + 1 ) 
*3-1 


03, 


03, 


(i-1) 

Stij 

*3+1 


CO. 


« \ 


iL 


03, 


w; 


(i-i) 

iL 

6 + 1 ) 


iL 


(4.4) 


(4.5) 



14 L) 

(L) 
03-■ 
*3- 

c 

,+) 

UJ ij+ 1 ' ' 

.XL) i 
^iL ) 




,,(!) 

' ' +/-1 

ft) 1 ’ 

,,(1) 

+■+1 

. ,(1) \ 

• ' +L 


: = 

6-i) 

w ii 

, ,6+i) 

w ii 

, ,6-i) 

• • +j-l 

, ,6+1) 

• • +3-1 

ftp 1 ’ 

4 +1> 

6-i) 

+3+1 

, ,6+i) 

+3+1 

, ,6-i) 

• • +L 

, ,6+i) 

• • +L 

(4.6) 


, ,+) 

\ +Vi 

, +) 


,+) 

+3+1 

.XL) , 

• • +L / 



(n) 


As far as formulae (4.5) and (4.6) are concerned it is worth stressing that columns con¬ 
taining terms co. 


are absent from the definition of the matrices and Hij. Therefore 


formula (4.3) expresses each single term £ij(Z) as a linear combination of Z and the 


partial derivative diZ. This particular type of relation will be the main ingredient for 
establishing the following theorem. 


where 


K, := 


and 


(co. 


(i) 

il 


VAi 


/ 


CO 


K, := 


(i) 

il 


, XL) 
\ W il 


CO. 


CO. 


CO. 


atisfies the system of equations 

det(Kj) v 
det(Kj) 


,(i) 5 (1) 

ii—1 "'i ^22+1 * * 

, , (1) \ 

• +L 1 

X L ) h (L) XL) 

ii—1 "'i ^22+1 * * 

, XL) 

■ U iL / 


(1) 

ii—1 


CO. 


C CO, 


( 1 ) 

22+1 


(+) 

ii—1 


C CO, 


(L) 

22+1 


OJ, 


,( 1 )\ 
iL 


CO, 


(L) 

i/L / 


(4.7) 


(4.8) 


(4.9) 


8 



























Proof. Substitute (4.3) into (3.4) with n = i. This procedure leaves us with the equation 
diZ = BiZ where 

L 


hfdetiWi) - 


B, := 


j=i 


(4,10) 


'I'rdi'j - y.y>,t.f//„i 


j'=i 


Next we use determinant expansion by minors to verify that the denominator of Bi 
corresponds to det(Kj) while the numerator equals det(Kj). □ 


Simplifications. The matrices Kj and Kj respectively defined by (4.8 ) an d (4.9 ) ar e 
given in terms of functions and h^ obtained through the relations (3.5) and (4.1). 


As a matter of fact we are interested only in the determinants of these matrices and 
that is still liable to simplifications. In this way it is convenient to rewrite the required 
determinants as 


det(Kj) = JJa 1 L (Ai -/rfc) JJ ^det(Yi) 

k =i k =i a ^ Ak Ai) 

k^i 

det(Kj) = a 1_L (Aj - /i k ) JJ b L _^ k —^ det(Yj) 

k =i fc=i a ^ Ak Ai > 

k^i 


(4,11) 


with matrices Y,- and Yj obtained respectively from Kj and Kj under the replacement 
uff H- . The functions are then given by 


*? 


u. 


(n) _ a (Aj — Aj) a(A n — Aj) -r-r . . TT ®(Afc — Aj) 

ij b{Xj - Aj) a(A n - Aj) W J ~~ ^ b{X k - A,) 


u. 


(n) 


®(A n Aj) 

c(A n Aj) 


| ^ b(\n fJ'k) 


k =1 


o(Aj A n ) 
b{\i A n ) 


2 L 


| ^ ®(A n Pk) J^Jf 


for j 7 ^ n, 


®(A k A n ) 


k =1 


fc=l ^ n ) 

k^i,n 


(4.12) 


The explicit evaluation of det(Yj) and det(Yj) for small values of L reveals that they are 
in fact multivariate polynomials and this seems to be the case for arbitrary values of L. 
Moreover, although we shall not present a proof here, det(Yj) is a multivariate polynomial 
of degree L — 1 in each variable separately while det(Yj) is of degree L — 2. The inspection 
of these polynomials for small values of L also indicates that they share no common 
zeroes for generic values of the variables A j. Consequently, the ratio det(Yj)/det(Yj) is a 
rational function preserving the polynomial degree of both numerator and denominator. 
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The partition function Z. Theorem [3] states that the partition function Z satisfies a 
simple system of linear first-order PDEs whose coefficients are expressed as determinants. 


Considering (4.7) and (4.11) we then have 


d 4 Z = 


det(Yj 
det (Y, 


■Z . 


(4.13) 


Eq. (4.13) has interesting consequences for functions Z G C 1 (C i ) and the solution can 
be obtained as follows. 


Corollary 1 . The solution of (4-13) with asymptotic behavior 


2~ L ( L ~i) c l u JJ A z 

3 = 1 


L—l 


as 


A j —y oo , 


(4.14) 


admits the family of representations 


Z = c 


L -1 


det(Yj 


(4.15) 


Proof. Restricting our attention to functions Z G C 1 (C L ) we can conclude that the 
residues of the RHS of ( 4. 13| ) at the zeroes of det(Yj) must vanish. In this way, Z needs 
to have the same zeroes as det(Yj) since det(Yj) does not share zeroes with det(Yj) for 
generic values of A j. We lack a proof for the latter statement but we have verified the 
validity of this property for several values of L. As previously mentioned, det(Yj) is a 
polynomial of degree L — 1 in each variable A j separately and consequently it has L — l 
zeroes with respect to a given variable A j. Therefore, Z needs to have at least L — l zeroes 
with respect to each variable A j. Now let V 3 := CfAf 1 ,..., A A 1 ,, AS X , ■ • •, Af 1 ] and define 


the set of zeroes Lf := {A G V 3 \ Z | A =A = 0} and Cj := {A G Vj j det(Y ?: )| A = A = 0}. 


Thus card(£j i ) = L — l while card (£f) > L — l. Next assume 3^ G and suppose 


C is an integration contour enclosing solely f. Under these assumptions Eq. (4.13) can 
be integrated as follows, 

M dAj /de tqq 

Z J J c det(Yi) 


'3 ' 


(4.16) 


The RHS of (4.16) vanishes while the LHS is different from zero if such element f exists. 


Hence we can conclude that Cf = Cj i by contradiction and consequently Z = ri det(Yj). 

By doing 

□ 


The constant r t is then fixed by the asymptotic behavior (4.14) derived in 
so we find r\ = c i_1 . 


31 


Remark 2. The RHS of (f.15) depends on an index i G {1, 2, ..., L} and it consequently 
comprises a total of L representations for the partition function Z. Although we have not 
presented a rigorous proof here, the explicit evaluation o/det(Yj) and det(Yj) fori j 
and several values of L confirms they result in the same partition function Z. Moreover, 
det(Yj) and det(Yj) do not seem to be related by any obvious determinant preserving 
transformation and this fact suggests they indeed consist of independent representations. 


10 

















5 Concluding remarks 


In this work we have studied the partition function of the six-vertex model with domain- 
wall boundaries from the perspective of partial differential equations. We have shown 
that a previously obtained functional equation describing the model’s partition function 
can be reduced to a system of first-order differential equations resembling the Knizhnik- 
Zamolodchikov equation. In fact, we have obtained a whole family of systems of equations 
which can be solved by elementary methods. The solution is then obtained as a single 
determinant whose partial homogeneous limit /jlj K > /j can be obtained trivially. Inter¬ 
estingly, the procedure described here yields a whole family of representations which do 
not seem to be related in any obvious way. 

The starting point of our analysis is the functional equation (2.1) which can also 
be solved in terms of a multiple contour integral as shown in |22,24|. However, the 


resolution of (2.1) along the lines described in 22 requires the use of additional properties 
of the desired solution, for example, they exhibit a particular polynomial structure. 
Furthermore, the implementation of the approach described in 22 depends strongly on 
the location of the so-called special zeroes which might not even exist for certain models. 
This seems to be the case for models based on non-compact symmetries such as the 
non-linear Schrodinger model 


38,39 


As for the method we are employing here to solve the system of PDEs (3.4), some 
remarks are also important. For instance, one of the arguments we have used for ob¬ 
taining (4.15) is that det(Y*) and det(Yj) are polynomials of particular degrees sharing 
no common zeroes. A rigorous proof of both arguments is missing but we have verified 
the validity of these assumptions for several lattice lengths. Nevertheless, under these 
assumptions we find that the full characterization of the desired partition function does 
not require much more input than what is already contained in our system of equations. 

The partition function of the six-vertex model with domain-wall boundary conditions 
has been already extensively studied in the literature and many of its properties are 
well known. For example, one of the key results in the literature is the determinant 
representation obtained in 14 which has played an important role for further studies of 
correlation functions of quantum integrable systems 39]. In spite of that, the structure 


of the functional equation (2.1) is shared by several other quantities as shown in 22 


35-37,40 and one can expect that our present analysis can be extended to those cases 


as well. More importantly, this analysis might be of particular relevance for the elliptic 
SOS model considered in 22 and 41 for which a single determinant representation for 


the model’s partition function is unknown to date. 
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